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T cell activation, initiated byT cell receptor (TCR) mediated recognition of pathogen-derived 
peptides presented by major histocompatibility complex class I or II molecules (pMHC), 
shows exquisite specificity and sensitivity, even though the TCR-pMHC binding interac- 
tion is of low affinity. Recent experimental work suggests that TCR pre-clustering may be a 
mechanism via whichT cells can achieve such high sensitivity. The unresolved stoichiometry 
of theTCR makesTCR-pMHC binding andTCR triggering, an open question. We formulate 
a mathematical model to characterize the pre-clustering of T cell receptors (TCRs) on the 
surface of T cells, motivated by the experimentally observed distribution of TCR clusters 
on the surface of naive and memory T cells. We extend a recently introduced stochas- 
tic criterion to compute the timescales of T cell responses, assuming that ligand-induced 
cross-linked TCR is the minimum signaling unit. We derive an approximate formula for the 
mean time to signal initiation. Our results show that pre-clustering reduces the mean activa- 
tion time. However, additional mechanisms favoring the existence of clusters are required 
to explain the difference between naive and memory T cell responses. We discuss the 
biological implications of our results, and both the compatibility and complementarity of 
our approach with other existing mathematical models. 

Keywords: T cell receptor, clustering, stochastic dynamics, signaling, naiveT cells, memoryT cells 



Edited by: 

Rob J. De Boer, Utrecht University, 
Netfierlands 

Reviewed by: 

Daniel Coombs, University of British 
Columbia, Canada 
Barbara Szomolay, University of 
Warwick, UK 

*Correspondence: 

Mario Castro, Universidad Pontificia 
Comiiias, C/ Alberto Aguilera 25, 
Madrid E28015, Spain 
e-mail: marioc@upcomillas. es; 
Hisse M. van Santen, Centre Biologia 
Molecular Severe Ochoa, Calle 
Nicolas Cabrera 7, Campus de 
Cantoblanco, Madrid 28049, Spain 
e-mail: hvansanten@cbm. csic. es 



1. INTRODUCTION 

A hallmark of the adaptive immune system is the ability of T cells, 
making use of the T cell receptors (TCRs) on their surface, to 
recognize a given agonist peptide-MHC ligand complex (pMHC) 
with high sensitivity (1). Some aspects of TCR-pMHC molecular 
interactions that are of current research interest are the frequency 
of encounters between T cells and the agonist pMHC, how cell- 
cell interactions determine the activation of lymphocytes (2), how 
early interactions change the state of the T cell receptor (3), what 
are the mechanisms of modulation of receptor-ligand interactions 
at cell-cell interfaces (4), and how protein organization in the cell 
membrane (for instance, protein islands or lipid rafts) affects the 
recognition process (3). Some recent experiments have explored 
the role of dimensionality on T cell activation and have highlighted 
the significance of the events taking place at the receptor level [see 
Refs. (1) and (6) for comprehensive reviews]. 

These open questions have been addressed with the use 
of mathematical modeling. Different theories can be classified 
according to the level of description (7). At the individual TCR- 
pMHC bond level, the kinetic proof-reading model (8) assumes 
that the TCR needs to undergo a series of consecutive (phos- 
phorylation) steps before being triggered. Also at the TCR level, 
the optimal dwell time model (9) reconciles the concurrence of 
different timescales, providing an optimal timescale between the 



very short times related to the off rate of TCR-pMHC binding, 
and the long times related to kinetic proof-reading mechanisms. 
The TCR occupancy model (10) considers the cell as a counting 
device in which multiple TCR-pMHC interactions are required 
to activate a T cell. In a similar fashion, the serial triggering 
model (II) proposed that the same pMHC can engage serially 
different TCRs. This model enriches the viewpoint of the TCR 
occupancy model, by giving greater relevance to the role of the 
pMHC itself Finally, the serial encounter model (12) and the con- 
finement time model (13) combine several of the ideas above and 
provide some appealing explanations by relaxing some restrictions 
in those models. 

While antigen presenting cells (APCs), such as dendritic cells 
or B cells, present 10^-10* times more self-pMHC than antigenic 
pMHC, self-pMHC ligands by themselves do not usually elicit a 
T cell response, even though their affinity for TCRa/3 is only 10 
times lower than the affinity of the antigenic pMHC (14). This 
illustrates how a small difference in affinity results in high speci- 
ficity, when there is only a few antigenic pMHC molecules in a 
background of self-pMHC ligands (15). 

The T cell signaling process begins with (extracellular) TCR- 
pMHC binding, followed by phosphorylation of the intracellular 
ITAM domains of the TCR-CD3 complex. When a TCR binds 
a pMHC molecule, the TCRa/J hetero-dimer binds the peptide. 
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while the CD4 or CDS co-receptor binds the MHC molecule. The 
binding of the co-receptor activates the tyrosine kinase Lck, which 
phosphorylates the ITAMs of the CD3 complex. ITAM phospho- 
rylation allows recruitment of intracellular signaling components 
that mediate downstream signaling events (16). 

It has recently been suggested that, contrary to what happens in 
TCR micro-clusters and the immunological synapse, clustering is 
not only induced by the ligand but by an avidity maturation mech- 
anism (or pre-clustering) (17), allowing the aggregation of chains 
of TCRs as long as 20 units (around 200 nm long), and referred to 
as nano-clusters (3, 18). Specifically, multimeric TCR-CD3 com- 
plexes are activated at low agonistic pMHC concentrations and 
monomeric TCRs remain unaffected at low ligand concentration. 
The TCR nano-clusters could enhance T cell sensitivity by the 
mechanisms proposed in the models of T cell activation (7), as 
their existence would reduce the time needed for two (or more) 
receptors to aggregate (by diffusion). This pre-cluster formation 
could be explained by three different mechanisms (3): 

• Multimeric complexes (or clusters) enhance the TCR avidity 
toward the ligand, which is expressed in clusters on the surface 
of APCs (19-21). At low ligand concentration, only multimeric 
TCR clusters are bound to ligand, as TCR monomers require 
higher ligand concentration. Monomeric TCRs might only be 
activated at high agonist doses. 

• Multimeric complexes allow the propagation of the activation 
signal from ligand-bound TCRayS to neighboring receptors in 
the same TCR multimer. 

• Linear arrays of multimeric TCR complexes help a single pMHC 
serially trigger several receptors (11). 

The existence of these nano-clusters does not exclude addi- 
tional mechanisms of T cell activation, as long as they involve the 
cooperation of receptors when they aggregate. Thus, while models 
such as kinetic proof-reading [and improvements as described in 
Ref. (22)] operate at the level of a single receptor, other models 
might be used in combination with the fact that the pre-cluster 
distribution of naive and memory T cells is different. 

Additionally, the fact that the TCR stoichiometry has not been 
resolved under physiological conditions, yet, makes it even more 
difficult to understand, at a molecular level, the dynamics of TCR 
pre-clustering (23). TCR pre-clustering could be an example of 
a more general mechanism of membrane-bound molecular pre- 
clustering, as clustering prior to cell-cell interaction has also been 
observed on the surface of APCs (19-21). It is worth mentioning 
that monomeric TCRs can still be activated at increasing lig- 
and concentrations, thus, conferring the T cell with a capacity 
to generate a dose-dependent response at very high pMHC doses, 
when multimeric TCR-CD3 complexes are already saturated (18). 
Such mechanisms have been previously described for chemo tactic 
bacteria, as a cellular mechanism to control sensitivity (24). 

Various mechanisms have already been suggested, at the pop- 
ulation, cellular or molecular level, to explain the capacity of T 
cells to respond, faster and more strongly, to a second antigenic 
encounter. However, the underlying mechanisms of the observed 
changes in the sensitivity of the T cell for pMHC ligand-mediated 



TCR stimulation (25) have not yet been clearly elucidated. Inter- 
estingly, the distribution of clusters in naive and memory T cells 
is different: memory T cells accommodate larger linear TCR clus- 
ters than naive ones. This could explain why memory T cells elicit 
more rapid responses than naive T cells (17) (see Figure 1 below). 

In this paper, we explore the consequences of TCR pre- 
clustering in signaling and in distinguishing naive from memory 
T cell responses. We present some experimentally obtained dis- 
tributions of TCR clusters for both types of cells (see Figure 1), 
and two complementary theoretical models: (i) a simple model 
of receptor oligomerization that describes cluster size distribu- 
tions, and (ii) a generalization of the stochastic T cell response 
criterion of Ref (26), to accommodate the hypothesis that the 
minimum signaling unit is composed of a TCR receptor cluster 
that is bound by the same cross-linked multivalent ligand. We 
find that this signaling unit is able to discriminate between ago- 
nist and antagonist pMHC ligands (with greater sensitivity than 
in the monomeric case), and to explain some of the advantages 
that higher cluster sizes can provide to memory T cells. The model 
also points at the need to invoke additional cooperativity mecha- 
nisms, to explain the experimentally observed role of clustering in 
T cell responses (27). Finally, this model of ligand-induced TCR 
cross-linking can be relevant in physiological conditions, accord- 
ing to the defective ribosomal products (DRiP) hypothesis (28,29), 
which provides a rapid source of peptide precursors to optimize 
immuno-surveillance of pathogens and tumors (30). 

2. MATHEMATICAL MODELING OF TCR PRE-CLUSTERING 
ANDT CELL ACTIV'OTini\! 

2.1. MODEL 1: T CELL RECEPTOR PRE-CLUSTERING 

The TCR-CD3 complex consists of the pMHC binding TCRa^S 
hetero-dimer, associated with the hetero-dimers CD3ye and 
CD336, and the homo-dimer CD3f f . Binding of a stimulating 
pMHC ligand by the extracellular domain of TCRa/J results in 
conformational changes in the intracellular part of the CD3e 
chain, and phosphorylation of the immuno-receptor tyrosine- 
based activation motifs (ITAMs) in the intracellular domains 
of the CD3y6, CD3i56, and CD3ff dimers, which in turn 
lead to initiation of downstream signaling cascades and T cell 
activation. 

It has long been recognized that the TCR-CD3 complex forms 
clusters upon ligand binding (31-36). More recently, it has been 
shown that in the absence of stimulating pMHC ligand, TCR- 
CD3 complexes are already expressed at the cell surface as a 
combination of monomeric and oligomeric TCR complexes or 
TCR nano-clusters (18). Electron microscopy (EM) analysis of 
immuno-gold-labeled human and murine T cells showed that 
these nano-clusters consist of up to 20 TCR-CD3 complexes. The 
exact stoichiometry of the nano-clusters has not been resolved yet. 

The integrity of TCR nano-clusters depends on cholesterol 
present at the cell surface membrane (18). The formation of the 
clusters depends, at least, on the trans-membrane region of the 
CD3ff homo-dimer (17), perhaps due to the capacity of f f 
dimers to form dimers of dimers (37). Other possible mecha- 
nisms of cluster formation rely on the capacity of the extracellular 
domain of TCRo; to dimerize (38). 
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FIGURE 1 I Distribution of TCRs at the surface of naive and memory 
T cells. Resting naive and memory CD8+ OT-1 T cells were labeled with the 

CD3e-specific mAb 2C11 and 10 nm gold-conjugated protein-A. Cell surface 
replicas of the labeled! cells were analyzed by transmission electron 
microscopy and the number and size of the observed gold clusters were 
recorded. (A)TEM image of surface replicas of a memory and a naive OT-1 T 
cell. The insets to the right show an enlargement of the boxed areas. 
(B) Quantification (mean± SD) of gold particles in clusters of the indicated 



sizes for resting naive T cells (gray bars, 7 cells, 9190 particles) and memoryJ 
cells (black bars, 5 cells, 3001 particles). The inset shows a detailed view of the 

distribution of clusters of three or more gold particles and statistical analysis 
(2-tailed Student's Most: •p<0.05, ••p< 0.01, and •••p< 0.001). All naive 
and memoryT cells had clusters with gold. However, whereas in naiveT cells 
the maximum gold cluster size shared by all cells was four, this was eight for 
memoryT cells. Also clusters bigger than twenty three particles were present 
in four out of five memoryT cells, and only two out of seven naiveT cells. 



www.frontiersin.org 



April 2014 I Volume 5 | Article 132 | 3 



Castro et al. 



TCR clustering andT cell response 



This body of experimental evidence allows us to conclude 
that multimeric TCR-CD3 complexes are co-expressed with TCR 
monomers on the surface of resting T cells. 

A simple model of aggregation of TCRa/J units is depicted in 
the left panel of Figure 2. Given a chain of length n (with n hetero- 
dimers linked), in a small time interval A t, with probability q+ A f, 
the chain increases to length n + 1 , and with probability q- At, the 
chain decreases to length n — I. Thus, by probability conservation, 
the probability to remain the same length « is 1 — (q+ + q_) At. 

Mathematically, the dynamics of the process can be described 
by a continuous time Markov chain (39) (or birth and death 
process, as we assume that polymerization takes place in unit 
steps). The state space is {1, 2, 3, . . ., n — 1, «, d + 1, . . .}, where 
the number denotes the number of TCRs in a cluster: 

q- q- q- q— q— q— 

1^2^3^---^M— l^ni=;M-|-l-- - . 

9+ <i+ q+ 9+ q+ 9+ 

The forward Kolmogorov equations for the probability of 
having a cluster of size n are given by (40) 



5„(f) 



dt 
dpiit) 



dt 



q+pn-i(t) + q-pn+i(t) - (q+ + q-)pn(t) , V« > 2, 

q-p2(t) - q+pi(t)- 



The stationary probability distribution is then given by 



lim p„(t) 



b"-'(l - b) 



(1 



b <l, ne {1,2,3, 



(1) 



with b = and n „ the probability (in thermodynamic equi- 
librium) to have a cluster of size n. When b < 1 (the number of 
clusters with a given size, n, decreases as n increases), and taking 
into account that peripheral T cells have around Nmux — 3 x 10* 
receptors, the latter expression can be further simplified to 



\l-b), b < I, n G {1,2,3, 



(2) 



2.2. MODEL 2: A BIVALENT MODEL FOR T CELL ACTIVATION 

The TCR-pMHC binding model introduced in Ref. (26) consid- 
ered monovalent pMHC ligands binding to TCR monomers on the 
surface of a T cell. Monovalent ligands have been reported to elicit 
a T cell response (41-43), but only when they are immobilized on 
a surface (which makes it difficult to assess whether they are truly 
monovalent or not) . Yet, multivalent receptor-ligand interactions 
are required to elicit T cell responses in both CD4+ and CD8+ T 
cells. In what follows, and supported by a body of experimental 
work (3, 24, 44), we adopt the hypothesis that the minimum acti- 
vating unit is a TCR-pMHC cross-linked dimeric complex (31, 
45-47). We make use of the binding model (Model 2) with pMHC 
dimers (ligands) and dimeric TCRs (receptors), described in the 
right panel of Figure 2. 

Gold-labeling experiments support the existence of nano- 
clusters with more than two TCRs, yet it can be shown (see Section 
5.2) that the key parameter of the mathematical model is the frac- 
tion of monomeric to multimeric TCR clusters. Thus, without loss 
of generality, we will assume that all TCR clusters are dimeric. 

The biochemical reactions encoded by the right panel of 
Figure 2 are as follows: 

• A (bivalent) ligand can bind a free receptor with monomeric 
binding reaction rates {kon and kog). Although not shown in 
the figure, we allow for a second ligand to bind the free receptor 
of the cluster. However, at low concentrations of ligands, this 
reaction can be safely neglected. 

• Cross-linking of a singly bound ligand follows with rates ki 
(forward reaction) and fc_2 (backward reaction). 

• If the complex formed by the ligand cross-linked to the dimeric 
TCR cluster lasts at least a time r, dwell time, we count that 
event. When we reach N such events, we will assume that a 
T cell response is initiated. The rationale behind this T cell 
response criterion follows the work of Ptilmer et al. (48), where 
the concepts of minimum dwell time and productive binding 
were introduced. This model combines aspects of the kinetic 
proof-reading (8) and the serial triggering models (7, 11). The 



Model 1 



n+1 



FIGURE 2 I Oligomerization and signaling models Left panel: 
oligomerization model not mediated by ligand (Model 1 ). We assume that 
receptors are able to diffuse and aggregate to an existing cluster. However, 
we exclude the possibility of clusters with size larger than one to diffuse. 
Clusters grow one monomeric unit at a time. Right panel: reactions included 
in the stochastic activation model (Model 2). Ligands in solution are able to 
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attach monovalently to any receptor in a cluster (top reaction). In addition, 
ligand-inducedTCR cross-linking can occur once a ligand is bound to a TCR 
in a given nano-cluster (bottom reaction). Following Ref. (26), once the 
bivalently bound ligand has been attached for a time z, we count that state 
as a signaling unit. After N of these units have been generated, the cell 
becomes activated. 
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minimum dwell time for a TCR-pMHC complex is the time the 
complex must remain bound in order to reach a level of ITAM 
phosphorylation, which will allow TCR triggering. Any binding, 
which persists for longer than the minimum dwell time is classi- 
fied as a productive binding [see Refs. (48) and (26) for further 
details] . 

• From an immunological perspective, the relevant parameter is 
the mean time to signal initiation, or MTSI (26). Namely, the 
MTSI is the average time needed for a T cell response according 
to the criterion that at least N TCR dimers should be bivalently 
engaged to a bivalent ligand (pMHC) for at least a time r. 

Here we assume that N is around 10-100. That is, 10-100 TCRs 
are required for signaling and = b x Nr is of the order of 10**, 
with Njj the total number of clusters on the T cell surface. This 
means, under the assumptions of Model 2, that at most, there can 
be N= 100 internalization events, as this is the number of trig- 
gered TCRs. Thus, in this approximation, the loss of TCR due to 
internalization after triggering can be safely neglected. Neverthe- 
less, internalization is an important step in early signaling, and 
a proper mechanistic model to justify the value of r will require 
internalization to be considered. This analysis is out of the scope 
of this article. 

We implement these reactions as a Markov process, and solve 
them numerically using the standard Gillespie algorithm (49), and 
with the parameters summarized in Table 1. We have made use of 
three different ligands: 4A, 4P, and 4N, which were also used in 
Ref. (26). For these ligands, that bind the same TCR with differ- 
ent affinities, a simple estimation of the number of cross-linking 
events required to elicit a T cell response is summarized in Table 2. 

There is some evidence that, under physiological conditions, the 
chance of two specific peptides being presented by two MHC mol- 
ecules in sufficient proximity and long enough to act as a dimer is 
very small (46). This will make ligand-induced TCR cross-linking 
a rare event. However, some recent experimental work on the dis- 
tribution of cognate pMHC molecules on the surface of APCs 
shows that both for MHC class I (virus infection models), and for 
MHC class II (antigen uptake via the endocytic route) clusters of 
cognate pMHC can be detected (19-21). 

We also note that ligand concentration is not the only fac- 
tor that depends on physiological conditions. According to the 
DRiP hypothesis (28,29), rapid viral antigen presentation is possi- 
ble because antigenic peptides originate from defective ribosomal 
products that have short half-lives. Although this phenomenon 
affects the time between viral challenge and antigen presentation, 
we assume it is independent of the subsequent signaling dynamics 
of T cell activation. 

3. RESULTS 

3.1 . DISTRIBUTION OF TCR CLUSTERS 

The mathematical model described in Section 2.1, or Model 1, 
allows us to obtain the value of h that best fits the experimental 
data. We have used a weighted (by the variance) minimum-square 
regression to fit the experimental distributions to equation (2). 
This kind of fit minimizes the value of x^- Thus, in Figure 3, we 
show the agreement between theory and experiment, with values: 
fcnaive = 0.32 and J'memoi-y = 0.55. The difference between bnaive 



Table 1 j Summary of the parameters used in the stochastic 
simulations. 



Parameter 


Value 


Comment 


Na 


6.023 X 10^3 


Avogadro's number 


Nr 


30,000 


Average number of TCRs perT cell (34) 


V 


50 /il 


Volume of the experiment 


Nc 


10^ cells 


Number of T cells in the experiment 


Vc 


V/Nc 


Average extracellular volume per cell 


k-2 


koff 


Cross4inking off rate 


ki 


koff(kd/k^'""") 


Cross4inking rate^ 


N 


10 


Minimum number of bound dimer-bivalent 






clusters to elicit aT cell response 


T 


1-4 s 


Dwell time 



For typical values of the dissociation rate, kj, we find that k^ is about 10-50 times 
k„ff. We have assumed k^2 = k^ following Fief (44). When not explicitly shown, 
we have used the same parameters as in Ref (26). 
"The cross-linking rate k2 is adapted from Ref. (44) for bivalent receptors. 

Table 2 | Estimated mean number of cross-linking events, 

N' ~ Ne""-^^ , required to elicit aT cell response (SP thymocytes). 



Ligand N' 





t{s) 


/V=10 


/V=100 


4P 


1 


3 


12 


{/Con = 153,691 M-i s-!) 


4 


3 


13 


(/toff = 0.01 69 s-i) 


8 


3 


14 


4A 


1 


7 


58 


= 157533 M-i s-i) 


4 


-10^ 


-10" 


{/Coff = 0.8664 s-i) 


8 


-10*5 


-10^ 


4N 


1 


-10^ 


~108 


(/Co„ = 149,385 M-i s""') 


4 


~1030 


~103i 


('^off=8.6e43 S"^ ) 


8 


~1080 


~106i 



and ^memory Can be explained by the existence of larger (or at least 
more localized) lipid rafts on the membrane of memory T cells (50, 
51). Thus, the rates q± could be the effective combination of two 
mechanisms: one related to the diffusion of receptors on the mem- 
brane, and the other related to the aggregation of the receptors at 
the molecular level. The presence of cholesterol on the membrane 
changes the diffusion coefficient of the TCR receptors, as receptor 
diffusion within the raft is inhibited due to protein anchorage (52) 
and, thus, stabilizes the formation of clusters (a larger value of b 
means that, once two receptors are embedded in the same lipid 
raft, it is more difficult for them to become separated from each 
other). 

A consequence of Model 1 is that, as the stationary probabilities 
need to sum up to one, the fraction of clusters of size larger than 
one is, precisely, b. This fraction is 72% higher for memory T cells 
than for naive T cells: bmemory/^naive = l-'72. 

3.2. MEAN TIME TO SIGNAL INITIATION 

In Figures 4A-D, we show how the stochastic criterion is able 
to provide a ligand hierarchy according to their potency. Namely, 
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FIGURE 3 I Comparison between the experimental distribution of 
clusters (see also Figure 1) and those from IVIodel 1 for (A) naiveT 
cells and (B) memoryT cells. The theoretical distribution has been 
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fitted to equation (2) using a weighted (by the variance) 
minimum-square regression. The fitted values are b„n„e = 0.32 and 
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FIGURE 4 I Dependence of the mean time to signal initiation 
(MTSI), T{N, t) to have N cross-linked ligand-receptor complexes 
bound for at least a dwell time t for different model parameters as 
shown in every panel. The results have been obtained by making use 
of a Gillespie algorithm, after averaging over 100 realizations for each 
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set of the parameters, summarized inTable 2 (a python code for the 
stochastic integration is available upon request). Units of time are 
seconds. All parameters are taken from Tables 1 and 2 except 
(A) W = 10, r =1 and /r2 = 10x <r„„; (B) W = 10, r =1 and /C2 = 40 x k„, 
(C) N= 10, T =4 and t-2 = 10x k„„; (D) A/=10, r =4 and *r2=40x k„„ 
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the most agonistic ligand, 4P, elicits a T cell response in times of 
the order of a few seconds in all cases. On the contrary, the most 
antagonistic ligand, 4N, takes extremely large times to do so (in 
practical terms, this means it does not elicit a T cell response). 
Thus, TCR clustering can enhance the potency of ligands, when 
compared to the monomeric case (26), as experimentally observed 
and theoretically shown. 

Following a similar approach to that of Ref (26), we can derive 
an approximate formula for the mean time to signal initiation 
(MTSI), T{N, t), for different ranges of ligand concentration, p. 
We write (see Figure 5A and Section 5.3 for further details): 



T{N,T) ~ 



T + 



T + 



Nexp(2<:-2T) 
4pNBK,n(k2/k-2) 



1/2 



at high concentration 
at intermediate 
concentration 

at low concentration 



(3) 



These three regimes correspond to different immunological 
scenarios. In the case of high concentration of ligand, ligand is 
in great excess, so that the required number of signaling units is 
reached, almost as soon as the first signaling unit is formed (time 
of order r). At low ligand concentration, the dynamics is lim- 
ited by the first binding event, as cross-linking occurs in a slower 
timescale. So, the MTSI has the same functional form as that for the 
monomeric case (26). Finally, for intermediate ligand concentra- 
tion, the competition between binding and cross-linking implies a 
more complicated mathematical relationship. Of greater relevance 
to the discussion is the nature of the ligand (with different fcon and 
koS rates), and the number of TCR clusters on the membrane of 
the T cell (encoded in the parameter Nb =bx Nr, with Nr, the 
average number of TCRs per T cell, see Table 1). 

An expression for the variance of the time to signal initiation 
(TSI) cannot be provided in a closed form [as done in Ref. (26)]. 
However, the fact that the variance decreases as the ligand concen- 
tration increases, suggests that the mathematical formula for the 
variance in the monovalent case can provide an upper bound to 
the present (dimeric) case. 



Using equation (3), we also can deduce the role of pre-clustering 
in the signaling time, or MTSI. As the number of bivalent clusters is 
b X Nr, the larger b is, the shorter the response time becomes. The 
model predicts that, for physiological conditions (not too high lig- 
and concentrations), the ratio of the MTSI for naive and memory 
T cells is inversely proportional to the ratio of their corresponding 
values of b. Namely, memory cells would respond up to 72% faster 
than naive ones (Figure 5B). 

4. DISCUSSION 

TCR triggering mechanisms are currently under debate [see, for 
example, Ref. (53) and (7) for recent reviews]. TCR clustering 
may be invoked as a description of the experimental results (27). 
The requirement for multivalent engagement of TCRs by pMHC 
ligands in CD4+ T cells has been widely shown (45, 47, 54, 55). 
The same requirement was shown in CD8+ T cells by Stone and 
Stern (56). 

In this paper, we have made use of the concept of mean time 
to signal initiation (MTSI or stochastic criterion) as a method 
to quantify the effect of TCR clustering on the timescales of T 
cell responses and, thus, to compare the behavior of naive and 
memory T cells. This criterion has also allowed us to compare 
the results in Section 3 for dimeric binding with those of Ref. 
(26) for monomeric binding. The introduction of the cross-linked 
ligand-receptor complex as the minimum signaling unit gives the 
response greater sensitivity to small differences in ligand affinity. 

A recent and novel feature of TCR immunology is the existence 
of TCR nano-clusters that are pre-formed, independently of ligand 
(3). This suggests that a simple stoichiometric clustering model 
(oligomerization of free TCRs diffusing on the T cell membrane) 
is enough to account for the distribution of TCR nano-clusters. 
In the case of naive T cells, Model 1 predicts an effective non- 
dimensional parameter, b= q-^^/q^, that allows us to explain the 
experimentally observed TCR cluster distributions. The presence 
of larger lipid rafts on the membrane of memory T cells might pro- 
vide support for the different values of b for naive and memory 
cells, bnaive and ??memory> respectively. It has recently been shown 
that receptor diffusion within the raft is inhibited due to protein 
anchorage (52). This reduction in the TCR diffusion coefficient 
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FIGURE 5 I (A) Comparison of the numerical solution of Model 2 (Gillespie algorithm with the parameters summarized inTable 2) and the approximate solution 
[equation (3)1 for ligand 4P and the same parameters as in Figure 4. (B) Comparison of the mean MTSI for naive (red dashed line) and memory (solid black line) 
T cells. 
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would increase the time required for the receptor to escape from 
the raft [in a similar fashion as other escape problems (37)]. This 
escape time is inversely proportional to the diffusion coefficient 
itself. A smaller TCR diffiisivity, as would be the case for memory 
T cells, will imply a larger residence time in the raft, which in turn 
will increase the probability of receptor aggregation in a given TCR 
cluster. A more detailed model of TCR diffusion and aggregation 
on the T cell membrane will be the subject of future work. 

Equation (3) shows the explicit dependence of the MTSI, T{N, 
t), on the parameter Ng, for given values of N and r. Nb is the 
average number of dimeric receptor clusters per T cell, so that 
NB=bx Nr, with Nr, the average number of TCRs per T cell 
(see Table 1). For large ligand concentration, the predicted T cell 
response time for memory and naive T cells is the same, and is 
equal to r. In the case of intermediate concentrations, the MTSI 
is proportional to . Finally, for low ligand concentration the 

MTSI is proportional to p This implies that, at low ligand con- 
centration, TCR pre-clustering alone can only account for at most 
72% of the reduction in the response time between memory and 
naive T cells. This behavior is illustrated in Figure 5A. This differ- 
ence is not so large as to be able to account for the observed higher 
responsiveness of memory T cells. Our results, thus, point to the 
need for additional mechanisms beyond TCR pre-clustering. 

A potential candidate to explain the large differences between 
memory and naive T cell responses is the conformational change 
of the CD3 complex (58). This conformational change is essen- 
tial to enable ITAM phosphorylation and, thus, the transfer of the 
TCR signal from the ecto-domain to the cytoplasmic tail of the 
TCR (58). Conformational changes in the CD3 complex occur as 
a result of the a/J hetero-dimer binding to pMHC. These con- 
formational changes allow the subunits of the CD3 complex (the 
ye and Se hetero-dimers and the f f homo-dimer) to become 
accessible to Lck, which can then phosphorylate their cytoplasmic 
domains at the ITAMs, leading to T cell signaling (59). In this 
way, the ligand-induced conformational change of the receptors 
can be propagated to all the receptors in the same cluster, so that 
larger clusters would benefit from this conformational change as a 
cascade [see, for example, Ref (60) and references therein]. Thus, 
differences in the distribution of cluster sizes could, indeed, explain 
the immunological differences between memory and naive cells. 

Other membrane receptors also exhibit pre-clustering and 
ligand-induced receptor cross-linking. For instance, in the case 
of the vascular endothelial growth factor receptor ( VEGFR), it has 
been shown (61) that there are two distinct pathways to receptor 
dimerization: (i) dynamic pre-dimerization (as the one described 
in Model 1), and (ii) ligand-induced receptor dimerization. The 
main conclusion in Ref (61) is that both mechanisms are almost 
indistinguishable at low ligand concentration. However, the first 
mechanism is more sensitive to changes in the binding affinity at 
large ligand concentration. Although the biological system stud- 
ied in Ref (61) is different from the T cell receptor considered 
here, their conclusions might be generalized as both receptors are 
tyrosine kinases. 

Bachmann et al. (62, 63) considered a model of diffusion and 
ligand-induced TCR clustering. Their model suggests that the exis- 
tence of large enough clusters greatly inhibits subsequent multimer 



diffusion, thus, reducing the relevance that this mechanism might 
have. This inhibition might be experimentally tested by exploiting 
the differences between naive (small and few clusters) and memory 
(large and many clusters). It will be interesting to make use of the 
models introduced in this paper to investigate the different roles of 
ligand binding and cellular activation (62), and TCR turnover (64). 

Finally, the existence of TCR pre-clusters [and the knowledge 
of their membrane distribution given by 7r„, equation (2)] can be 
considered in the kinetic-segregation model (65). In this model, 
diffusion out of close-contact zones would be inhibited by the exis- 
tence of nano-clusters, thus, enhancing the number of triggered 
receptors. In a similar way, consecutive receptor phosphorylation 
events (66) in TCR nano-clusters would also amplify receptor 
signaling. 

5. MATERIALS AND METHODS 

5.1. EXPERIMENTS 

Naive CD8+ OT-1 T cells, which recognize an ovalbumin-derived 
peptide presented by the MHC class I molecule H-2K'', were 
isolated from superficial and mesenteric lymph nodes of OT-1 
TCR transgenic mice (67), via depletion of CD19+ B cells, CD4+ 
helper T cells and CDllb+ macrophages, using antibodies and 
Dynal magnetic beads (Invitrogen). Memory OT-1 T cells were 
generated by adoptively transferring 10* naive OT-1 T cells into 
congenic C57BL/6 LyS.l Pep3b mice, which were simultaneously 
immunized with 10^ PFU MVA-OVA (68). After 6 months, rest- 
ing memory OT-1 T cells were isolated from the spleen and 
lymph nodes of these mice by antibody-mediated depletion of 
macrophages, B cells, and CD4+ T cells, followed by separation of 
the OT-1 memory T cells from host-derived Ly5.1+ CD8+ T cells 
via fluorescence-activated cell sorting, using a LyS.l -specific anti- 
body. Labeling of cells with the CD3e-specific antibody 2C11 and 
10 nm gold-conjugated protein-A, replica generation and analysis 
were performed as previously described (17). 

5.2. MODELS OF SIGNALING WITH DIMERIC AND TRIMERIC 
RECEPTOR CLUSTERS 

In Section 2, we introduced a model in which ligands are biva- 
lent and receptor clusters are dimeric (that is, composed of two 
monomeric TCRs). This is, of course, a first approximation that 
neglects the distribution of cluster sizes experimentally observed. 
Yet, the results of our mathematical study only change in a quanti- 
tative way, but not qualitatively, when we include TCR clusters of 
larger sizes. In this Section, we illustrate this by considering a sys- 
tem in which clusters of size 1, 2, and 3 coexist and the ligands are 
bivalent. Table 3 provides the notation introduced to describe the 
molecular species considered in the model, as well as a graphical 
representation. 

At large initial ligand concentration, under the stochastic crite- 
rion, the MTSI tends to r. On the other hand, at low initial ligand 
concentration, the number of receptors, compared to the number 
of ligands, is so large that we can neglect molecular species X4,y4, 
76, and yy, which involve more than one bivalent ligand. This has 
also been confirmed experimentally. Given our stochastic T cell 
response criterion, in this case, the signaling units correspond to 
molecular species xs, ys, and yy. Molecular species zi and zj do 
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Table 3 | Summary of variables for a model in which clusters of size 1, 
2, and 3 coexist. 



to 73 is that of the initial ratio of free TCR dimers to free TCR 
trimers, namely, 



Variable Description 



^1 

^2 
Z3 

Xl 
X3 

xa 

^6 

Kl 
Y2 

K3 

Yi 
Yb 

Y& 

K7 



Free monomeric receptor 
Free ligand (dimer) 

Ligand-bound to a monomeric receptor 
Free dimeric cluster 

Same as Z2 (defined for convenience of 
notation) 

Ligand singly bound to a dimeric cluster 

Two ligands bound to a dimeric cluster 
Cross-linked ligand in a dimeric cluster 

Free trimeric receptor 

Same as Z2 (defined for convenience of 
notation) 

Ligand singly bound to a trimeric cluster 

Two ligands bound to a trimeric cluster 
Cross-linked ligand in a trimeric cluster 

Three ligands bound to a trimeric cluster 

One ligand singly bound to a trimeric 
cluster and another cross-linked 



Molecular 
representation 



I I I 
I 

I I I I 

I 

1^1 I 



iSllii 



All the variables correspond to the total number of molecular species (not 
concentrations). Hence, all the rates In the mathematical model have units ofs^'. 

not contribute to the T cell response and will be neglected in what 
follows. Thus, we only need to consider the dynamics of dimeric 
and trimeric T cell receptor clusters. 

We introduce the total number of signaling units, Ssit) = 
^sit) + ysit)- The set of ordinary differential equations for the 
model is given by: 



X3 



^2 



73 - & ^3 , 



(4) 



where we have made use of equation (2) to conclude ^ = b. Thus, 
the total number of signaling units, Ssit), obeys the following 
differential equation 



Ss = X5 + 75 = kl (1 + 2b) X3 - 2 fc_2 S5 



(5) 



Finally, in the low ligand concentration limit as above, let us 
introduce S3 = X3 + 73. It is easy to show that equation (5) 
reduces to 



,l + 2fo 

S5 = k2 ——r S3 -2 k-2 S5 , 

1 + 0 



(6) 



which is identical to the differential equation for xs above, but with 
S5,3 replaced by , respectively. This means that, except for a pre- 
factor ^jt^ [which, for b e (0, 1), is between 1 and 3/2], the study 
of dimeric and trimeric clusters is reduced to the dimeric case. 



5.3. A SIMPLE FORMULA FOR THE MTSI 

The basic idea behind the stochastic criterion is to count the cumu- 
lative number of events that may contribute to signaling (26). Here, 
we calculate the mean number of cross-linking events up to time 
t, C{t), as the integral. 



C(f) 



Jo 



h / Xiis) ds . 
lo 



(7) 



It is possible to obtain an expression for x^it) with the approx- 
imation that the product xiit)x2it) is constant, so that the pair 
of equations for x^it) and xs(t) can be solved exactly. This 
yields (69): 



C{t) 



^(Ai+2fc_2)(e'"' 
^1 



+ ^(A2 + 2/:_2)(e''' 



- 1) 
1) + ait 



(8) 



xi = -4fc+xiX2 + kosXi, 

X2 = -Ak+XiX2 + /CoffX3, 

X3 = 4k^XiX2 — kgf(X3 — k2X3 + 2fc_2X5, 

X5 = /C2X3 - 2/c_2X5, 

71 = -6fc+7i72 + /Coff73, 

72 = -6L1-7172 + kosyi, 

73 = 6fc+7i72 - /Coff73 - 2k2y3 + 2fc_275> 



where 



2k 



■273 ■ 



2k^ 



275> 



where fc+ = kon/(VNA), V is the volume of the experiment and 
Na is Avogadro's number. 

Given the symmetry of the problem, and in the limit of low 
initial ligand concentration, we will assume that the ratio of X3 



-4X2konPNB 



C2 



fl2 



(>^2 - Ai)(4/ConPfa + 2fcoff A:_2 + 8fconP^2) 

4:XikonpNB 

{>^2 - Xi){AkonPk2 + 2/Coff /c_2 + 8konPk-2) ' 

^ (-4fconP - fcoff - k2- 2/c_2 

±[(4fconP + fcoff -h- 2k^2f + '^KiM^'^) 

8/c_2AonPNB 

4fconP^2 + 2kosk-2 + SkonPk-2 ' 

4/C2/ConpNB 

^K>nPh + 2fcoff^-2 + 8/ConpA:_2 ' 
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and Nb is the number of dimeric receptors. The MTSI is 
then given by the solution of the equation C{T{N, t)-t) = 
Nexp(2fc_2T)- 

The expressions in equation (3) are obtained from equa- 
tion (8) in the appropriate regimes. At low ligand concentra- 
tion, C(f) is simply proportional to time: C(f) — k2a\t, so 
that C(r-r) = A:2ai(r-T) = Nexp(2A:_2r). When Ai,2T « 1, 
on the other hand, the first non-zero term in a Taylor expansion of 
C(f) in time is quadratic: C(f) oc t^. This provides the exponent 
1/2 in the second line of equation (3). 
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